Science Bulletin

Neuroanatomical tracing technology is fundamental for unraveling the complex network of brain connectome. Tracing tools that could spread between neurons are urgently needed, especially the rigorous trans-monosynaptic anterograde tracer is still lacking. HSV1 strain H129 was proved to be an anterograde tracer and has been used to trace neuronal networks in several reports. However, H129 has a serious defect that it was demonstrated to infect neurons via axon terminals. Thus, when using H129 to dissect output neural circuit, its terminal take up capacity should be carefully considered. Here, we report a recombinant H129 that carrying the anti-Her2 scFv in glycoprotein D to target genetically defined neurons. With the usage of helper virus complementarily expressing Her2 and gD, we can realize the elucidation of direct projection regions of either a given brain nucleus or a specific neuron type. The retargeted H129 system complements the current neural circuit tracer arsenal, which provides a rigorous and practical anterograde trans-monosynaptic tool.

Hypertrophic cardiomyopathy (HCM) is the most common inherited heart disease. We performed a comprehensive analysis to construct the correlation of m6A and immune in HCM. Two HCM datasets (GSE141910 and GSE160997) and m6A-related regulators were obtained from GEO and published articles, respectively. Differentially expressed m6A-related regulators were obtained. Random forest model and nomogram were conducted to assess the risk of HCM, and finally, the m6A subtype was constructed. Functional enrichment analysis was conducted. Protein-protein interaction network of differentially expressed genes between m6A subtypes was performed. Furthermore, we constructed the Hubgene-chemical network, Hubgene-microRNA network, and Hubgene-transcription factor network of the top 10 hubgenes. Additionally, the immune subtype and hubgene subtype were constructed. PCR was performed to validate the m6A-related regulators. We obtained 20 m6A-related regulators in HCM. Among them, 8 m6A-related regulators differentially expressed (YTHDC1, HNRNPC, and FMR1 were up-regulated while YTHDC2, FTO, WTAP, IGF2BP2, and IGF2BP3 were down-regulated). FTO, FMR1, IGF2BP3, YTHDC1, and IGF2BP2 were the top 5 important m6A-related regulators and were used to conduct the nomogram. We obtained 329 differentially expressed genes in m6A subtype and these genes enriched HCM-related processes and pathways. Furthermore, we constructed the Hubgene-chemical network, Hubgene-microRNA network, and Hubgene-transcription factor network of the top 10 hubgenes (NFKBIA, NFKB1, PSMA3, PSMC4, PSMA2, PSMA4, PSMD7, PSMD10, PSMD8, and PSMA6). And then we constructed an immune subtype based on the immune cell infiltration levels and hubgene subtype based on the expression of the top 10 hubgenes. Finally, we verified the main results through experiments. In conclusion, we built a nomogram and identified 8 m6A-related regulators and 10 hubgenes, which were prominently associated with HCM. We found that m6A and the immune system may play a crucial role in the HCM. Accordingly, those genes and pathways might become therapeutic targets with clinical usefulness in the future.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) cause the most serious pandemics of Coronavirus Disease 2019 (COVID-19), which threatens human health and public safety. SARS-CoV-2 spike (S) protein uses angiotensin-converting enzyme 2 (ACE2) as recognized receptor for its entry into host cell that contributes to the infection of SARS-CoV-2 to hosts. Using computational modeling approach, this study resolved the evolutionary pattern of bonding affinity of ACE2 in 247 jawed vertebrates to the spike (S) protein of SARS-CoV-2. First, high-or-low binding affinity phenotype divergence of ACE2 to the S protein of SARS-CoV-2 has appeared in two ancient species of jawed vertebrates, Scyliorhinus torazame (low-affinity, Chondrichthyes) and Latimeria chalumnae (high-affinity, Coelacanthimorpha). Second, multiple independent affinity divergence events recur in fishes, amphibians-reptiles, birds, and mammals. Third, high affinity phenotypes go up in mammals, possibly implying the rapid expansion of mammals might accelerate the evolution of coronaviruses. Fourth, we found natural mutations at eight amino acid sites of ACE2 can determine most of phenotype divergences of bonding affinity in 247 vertebrates and resolved their related structural basis. Moreover, we also identified high-affinity or low-affinity-associated concomitant mutation group.The group linked to extremely high affinity may provide novel potentials for the development of human recombinant soluble ACE2 (hrsACE2) in treating patients with COVID-19 or for constructing genetically modified SARS-CoV-2 infection models promoting vaccines studies. These findings would offer potential benefits for the treatment and prevention of SARS-CoV-2.

The outbreak of COVID-19 across the world has posed unprecedented and global challenges on multiple fronts. Most of the vaccine and drug development has focused on the spike proteins and viral RNA-polymerases. Using bioinformatics and structural modeling approach, we modeled the structure of the envelope (E)-protein of novel SARS-CoV-2. The E-protein of this virus shares sequence similarity with that of SARS-CoV-1, and is highly conserved in the N-terminal regions. Incidentally, compared to spike proteins, E proteins demonstrate lower disparity and mutability among the isolated sequences. Using homology modeling, we found that the most favorable structure could function as a gated proton channel. Combining pocket estimation and docking with water, we determined that GLU 8 and ASN 15 in the N-terminal region were in close proximity to form H-bonds. Additionally, two distinct “core” structures were visible, the hydrophobic core and the central core, which may regulate the opening/closing of the channel. We propose this as a mechanism of viral proton channeling activity which may play a critical role in viral infection. In addition, it provides a structural basis and additional avenues for generating therapeutic interventions against the virus.One Sentence Summary Structural modeling of the novel coronavirus envelope proteins (E-proteins) demonstrating its possible proton channeling activity.Competing Interest StatementThe authors have declared no competing interest.View Full Text

Unsupervised AI (artificial intelligence) can obtain novel knowledge from big data without particular models or prior knowledge and is highly desirable for unveiling hidden features in big data. SARS-CoV-2 poses a serious threat to public health and one important issue in characterizing this fast-evolving virus is to elucidate various aspects of their genome sequence changes. We previously established unsupervised AI, a BLSOM (batch-learning SOM), which can analyze five million genomic sequences simultaneously. The present study applied the BLSOM to the oligonucleotide compositions of forty thousand SARS-CoV-2 genomes. While only the oligonucleotide composition was given, the obtained clusters of genomes corresponded primarily to known main clades and internal divisions in the main clades. Since the BLSOM is explainable AI, it reveals which features of the oligonucleotide composition are responsible for clade clustering. The BLSOM has powerful image display capabilities and enables efficient knowledge discovery about viral evolutionary processes.

The outbreak of COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, rapidly spread to create a global pandemic and has continued to spread across hosts from humans to animals, transmitting particularly effectively in mink. How SARS-CoV-2 evolves in animals and humans and the differences in the separate evolutionary processes remain unclear. We analyzed the composition and codon usage bias of SARS-CoV-2 in infected humans and animals. Compared with other animals, SARS-CoV-2 in mink had the most substitutions. The substitutions of cytidine in SARS-CoV-2 in mink account for nearly 50% of the substitutions, while those in other animals represent only 30% of the substitutions. The incidence of adenine transversion in SARS-CoV-2 in other animals is threefold higher than that in mink-CoV (the SARS-CoV-2 virus in mink). A synonymous codon usage analysis showed that SARS-CoV-2 is optimized to adapt in the animals in which it is currently reported, and all the animals showed decreased adaptability relative to that of humans, except for mink. A binding affinity analysis indicated that the spike protein of the SARS-CoV-2 variant in mink showed a greater preference for binding with the mink receptor ACE2 than with the human receptor, especially as the mutation Y453F and F486L in mink-CoV lead to improvement of binding affinity for mink receptor. Our study focuses on the divergence of SARS-CoV-2 genome composition and codon usage in humans and animals, indicating possible natural selection and current host adaptation.

Mental processing delineates the functions of human mind encompassing a wide range of motor, sensory, emotional, and cognitive processes, each of which is underlain by the neuroanatomical substrates. Identifying accurate representation of functional neuroanatomy substrates of mental processing could inform understanding of its neural mechanism. The challenge is that it is unclear whether a specific mental process possesses a functional neuroanatomy fingerprint, i.e., a unique and reliable pattern of functional neuroanatomy that underlies the mental process. To address this question, we utilized a multi-task deep learning model to disentangle the functional neuroanatomy fingerprint of seven different and representative mental processes including Emotion, Gambling, Language, Motor, Relational, Social, and Working Memory. Results based on the functional magnetic resonance imaging data of two independent cohorts of 1235 subjects from the US and China consistently show that each of the seven mental processes possessed a functional neuroanatomy fingerprint, which is represented by a unique set of functional activity weights of whole-brain regions characterizing the degree of each region involved in the mental process. The functional neuroanatomy fingerprint of a specific mental process exhibits high discrimination ability (93% classification accuracy and AUC of 0.99) with those of the other mental processes, and is robust across different datasets and using different brain atlases. This study provides a solid functional neuroanatomy foundation for investigating the neural mechanism of mental processing. One-Sentence SummaryThere exists a functional neuroanatomy fingerprint to underlie the mental process.

The SARS-CoV-2 infected cases and the caused mortalities have been surging since the COVID-19 pandemic. Viral mutations emerge during the virus circulating in the population, which is shaping the viral infectivity and pathogenicity. Here we extensively analyzed 6698 SARS-CoV-2 whole genome sequences with specific sample collection dates in NCBI database. We found that four mutations, i.e., 5UTR_c-241-t, NSP3_c-3037-t, NSP12_c-14408-t, and S_a-23403-g, became the dominant variants and each of them represented nearly 100% of all virus sequences since the middle May, 2020. Notably, we found that co-occurrence rates of three significant multi-site co-mutational patterns, i.e., (i) S_a-23403-g, NSP12_c-14408-t, 5UTR_c-241-t, NSP3_c-3037-t, and ORF3a_c-25563-t; (ii) ORF8_t-28144-c, NSP4_c-8782-t, NSP14_c-18060-t, NSP13_a-17858-g, and NSP13_c-17747-t; and (iii) N_g-28881-a, N_g-28882-a, and N_g-28883-c, reached 66%, 90%, and nearly 100% of recent sequences, respectively. Moreover, we found significant decrease of CpG dinucleotide at positions 241(c)-242(g) in the 5UTR during the evolution, which was verified as a potential target of human zinc finger antiviral protein (ZAP). The four dominant mutations, three significant multi-site co-mutations, and the potential escape mutation of ZAP-target in 5UTR region contribute to the rapid evolution of SARS-CoV-2 virus in the population, thus shaping the viral infectivity and pathogenicity. This study provides valuable clues and frameworks to dissect the viral replication and virus-host interactions for designing effective therapeutics. One Sentence SummaryFour dominant mutations, three significant multi-site co-mutations, and 5UTR CpG escape contribute to the rapid evolution of SARS-CoV-2 virus.

Preventing immune escape of SARS-CoV-2 variants is crucial in vaccine development to ensure broad protection against the virus. Conformational epitopes beyond the RBD region are vital components of the spike protein but have received limited attention in the development of broadly protective SARS-CoV-2 vaccines. In this study, we used a DNA prime-protein boost regimen to evaluate the broad cross-neutralization potential of immune response targeting conformational non-RBD region against SARS-CoV-2 viruses in mice. Mice with enhanced antibody responses targeting conformational non-RBD region show better performance in cross-neutralization against the Wuhan-01, Delta, and Omicron subvariants. Via analyzing the distribution of conformational epitopes, and quantifying epitope-specific binding antibodies, we verified a positive correlation between the proportion of binding antibodies against the N-terminal domain (NTD) supersite (a conformational non-RBD epitope) and SARS-CoV-2 neutralization potency. The current work highlights the importance of high ratio of conformational non-RBD-specific binding antibodies in mediating viral cross-neutralization and provides new insight into overcoming the immune escape of SARS-CoV-2 variants.

Pandemic coronavirus disease 2019 (COVID-19) is caused by the emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), for which there are no efficacious vaccines or therapeutics that are urgently needed. We expressed three versions of spike (S) proteins--receptor binding domain (RBD), S1 subunit and S ectodomain--in insect cells. RBD appears monomer in solutions, whereas S1 and S associate into homotrimer with substantial glycosylation. The three proteins confer excellent antigenicity with six convalescent COVID-19 patient sera. Cryo-electron microscopy (cryo-EM) analyses indicate that the SARS-CoV-2 S trimer dominate in a unique conformation distinguished from the classic prefusion conformation of coronaviruses by the upper S1 region at lower position ~15 [A] proximal to viral membrane. Such conformation is proposed as an early prefusion state for the SARS-CoV-2 spike that may broaden the knowledge of coronavirus and facilitate vaccine development.

For adult diffuse gliomas (ADGs), most grading can be achieved through molecular subtyping, retaining only two key histopathological features for high-grade glioma (HGG): necrosis (NEC) and microvascular proliferation (MVP). We developed a deep learning (DL) framework to automatically identify and characterize these features. We trained patch-level models to detect and quantify NEC and MVP using a dataset that employed active learning, incorporating patches from 621 whole-slide images (WSIs) from the Chinese Glioma Genome Atlas (CGGA). Utilizing trained patch-level models, we effectively integrated the predicted outcomes and positions of individual patches within WSIs from The Cancer Genome Atlas (TCGA) cohort to form datasets. Subsequently, we introduced a patient-level model, named PLNet (Probability Localization Network), which was trained on these datasets to facilitate patient diagnosis. We also explored the subtypes of NEC and MVP based on the features extracted from patch-level models with clustering process applied on all positive patches. The patient-level models demonstrated exceptional performance, achieving an AUC of 0.9968, 0.9995 and AUPRC of 0.9788, 0.9860 for NEC and MVP, respectively. Compared to pathological reports, our patient-level models achieved the accuracy of 88.05% for NEC and 90.20% for MVP, along with a sensitivity of 73.68% and 77%. When sensitivity was set at 80%, the accuracy for NEC reached 79.28% and for MVP reached 77.55%. DL models enabled more efficient and accurate histopathological image analysis which will aid traditional glioma diagnosis. Clustering-based analyses utilizing features extracted from patch-level models could further investigate the subtypes of NEC and MVP.

Neurotropic virus tracers, particularly those with low toxicity and high efficient tracing, are powerful tools for structural and functional dissections of neural circuits. The retrograde trans-mono-synaptic technology based on rabies virus CVS-N2c strain has reduced cytotoxicity and enhanced efficiency, attains long-term gene manipulation for functional studies, but suffers from difficult preparation and low yield. To overcome these shortcomings, an improved production system was established for rapid rescue and preparation of CVS-N2c-{Delta}G virus, CVS-N2c-{Delta}G with the same titer as SAD-B19-{Delta}G can be prepared within a short time. Meanwhile, we found that N2cG coated CVS-N2c-{Delta}G allows efficient retrograde access to projection neurons, and further expand its application in VTA/SNc to DLS pathway that unaddressed by rAAV9-Retro, and the efficiency is 6 folds higher than that of rAAV9-Retro. Then the trans-synaptic efficiency of CVS-N2c-{Delta}G virus was evaluated. Results showed that the trans-mono-synaptic efficiency of oG-mediated CVS-N2c-{Delta}G was 2-3 folds higher than that of oG-mediated SAD-B19-{Delta}G, but there was no difference between oG-mediated and N2cG-mediated CVS-N2c-{Delta}G system. In addition, codon modified N2cG (optiG) did not increase the efficiency of CVS-N2c-{Delta}G tracing. Finally, we found that the CVS-N2c-{Delta}G produced by the improved method can be used for monitoring neural activity of projection neurons, and the time window can be maintained for 3 weeks, and it can also express sufficient recombinases for efficient transgene recombination. That is, the virus produced by the improved production system does not affect its own function, paving the way for its further optimization, popularization and application in structural and functional studies of neural circuits.

Serine/Arginine-Rich Splicing Factor 7 (SRSF7), which is previously recognized as a splicing factor, has been revealed to play oncogenic roles in multiple cancers. However, the mechanisms underlying its oncogenic roles have not been well addressed. Here, based on N6-methyladenosine (m6A) co-methylation network analysis across diverse cell lines, we found SRSF7 positively correlated with glioblastoma cell-specific m6A methylation. We then proved SRSF7 is a novel m6A regulator that specifically facilitates the m6A methylation near its binding sites on the mRNAs involved in cell proliferation and migration through recruiting methyltransferase complex. Moreover, SRSF7 promotes the proliferation and migration of glioblastoma cells largely dependent on the m6A methyltransferase. The two single-nucleotide m6A sites on PBK are regulated by SRSF7 and partially mediate the effects of SRSF7 on glioblastoma cells through recognition by IGF2BP2. Together, our discovery revealed a novel role of SRSF7 in regulating m6A and timely confirmed the existence and functional importance of RNA binding protein (RBP) mediated specific regulation of m6A.

Three-dimensional (3D) histopathology involves the microscopic examination of a specimen, which plays a vital role in studying tissues 3D structures and the signs of diseases. However, acquiring high-quality histological images of a whole organ is extremely time-consuming (e.g., several weeks) and laborious, as the organ has to be sectioned into hundreds or thousands of slices for imaging. Besides, the acquired images are required to undergo a complicated image registration process for 3D reconstruction. Here, by incorporating a recently developed vibratome-assisted block-face imaging technique with deep learning, we developed a pipeline termed HistoTRUST that can rapidly and automatically generate subcellular whole organs virtual hematoxylin and eosin (H&E) stained histological images which can be reconstructed into 3D by simple image stacking (i.e., without registration). The performance and robustness of HistoTRUST have been successfully validated by imaging all vital mouse organs (brain, liver, kidney, heart, lung, and spleen) within 1-3 days depending on the size. The generated 3D dataset has the same color tune as the traditional H&E stained histological images. Therefore, the virtual H&E stained images can be directly analyzed by pathologists. HistoTRUST has a high potential to serve as a new standard in providing 3D histology for research or clinical applications.

Host-virus protein-protein interaction is the key component of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lifecycle. We conducted a comprehensive interactome study between the virus and host cells using tandem affinity purification and proximity labeling strategies and identified 437 human proteins as the high-confidence interacting proteins. Functional characterization and further validation of these interactions elucidated how distinct SARS-CoV-2 viral proteins participate in its lifecycle, and discovered potential drug targets to the treatment of COVID-19. The interactomes of two key SARS-CoV-2 encoded viral proteins, NSP1 and N protein, were compared with the interactomes of their counterparts in other human coronaviruses. These comparisons not only revealed common host pathways these viruses manipulate for their survival, but also showed divergent protein-protein interactions that may explain differences in disease pathology. This comprehensive interactome of coronavirus disease-2019 provides valuable resources for understanding and treating this disease.

Formalin-fixed and paraffin-embedded (FFPE) tissue archives are the largest repository of clinically annotated specimens, and FFPE-compatible single cell gene expression workflow had been developed and applied recently. However, for tissues where cells are hard to dissociate or brains with complex neuronal cells, nuclear transcriptomic profiling are desirable. Moreover, the effects of standard pathological practice on the transcriptome of samples obtained from such archived specimens was also largely anecdotal. Here, we performed RNA-seq of nuclei from hippocampal of mice that underwent freezing, paraformaldehyde (PFA) fixation, and paraffin embedding. Then, we comprehensively evaluated the parameters affecting mRNA quality, transcription patterns, functional level and cell states of nuclei, including PFA fixation time and storage time of FFPE tissues. The results showed that the transcriptome signatures of nuclei isolated from fresh PFA-fixed and fresh FFPE tissues were more similar to matched frozen samples. By contrast, the brain fixed for more than 3 days had prominent impacts on the sequencing data, such as the numbers and biotypes of gene, GC content and ratio of reads interval. Commensurately, prolonged fixation time will result in more differentially expressed genes, especially those enriched in spliceosome and synaptic related pathways, affecting the analysis of gene splicing and neuron cells. MuSiC deconvolution results revealed that PFA infiltrating brains for 3 days will destroy the real cell states, and the proportion of neuron, endothelial and oligodendrocytes diminished while that of microglia was reversed. Yet the effect of storage time on cell composition was more neglectable for FFPE samples. In addition, oligodendrocyte precursor cells were most affected in all fixed samples, and their destruction was independent of fixation time and preservation time. The comprehensive results highlighted that fixation time had much more influences on the nuclear transcriptomic profiles than FFPE retention time, and the cliff-like effects appeared to occur over a fixed period of 1-3 days, with no more differences from additional fixation durations.

The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) becomes a tremendous threat to global health. Although vaccines against the virus are under development, the antigen epitopes on the virus and their immunogenicity are poorly understood. Here, we simulated the three-dimensional structures of SARS-CoV-2 proteins with high performance computer, predicted the B cell epitopes on spike (S), envelope (E), membrane (M), and nucleocapsid (N) proteins of SARS-CoV-2 using structure-based approaches, and then validated the epitope immunogenicity by immunizing mice. Almost all 33 predicted epitopes effectively induced antibody production, six of which were immunodominant epitopes in patients identified via the binding of epitopes with the sera from domestic and imported COVID-19 patients, and 23 were conserved within SARS-CoV-2, SARS-CoV and bat coronavirus RaTG13. We also found that the immunodominant epitopes of domestic SARS-CoV-2 were different from that of the imported, which may be caused by the mutations on S (G614D) and N proteins. Importantly, we validated that eight epitopes on S protein elicited neutralizing antibodies that blocked the cell entry of both D614 and G614 pseudo-virus of SARS-CoV-2, three and nine epitopes induced D614 or G614 neutralizing antibodies, respectively. Our present study shed light on the immunodominance, neutralization, and conserved epitopes on SARS-CoV-2 which are potently used for the diagnosis, virus classification and the vaccine design tackling inefficiency, virus mutation and different species of coronaviruses.

AimsTo analyze the role of lncRNA PVT1 in cerebral ischemia-reperfusion (I/R) injury induced by deep hypothermia low flow (DHLF). Methods and resultsA total of 72 lncRNAs were differentially expressed in the brain tissue of DHLF mice. PVT1 expression was significantly downregulated in DHLF mouse brain tissue, preoperative and postoperative blood samples from children undergoing DHLF extracorporeal circulation, and hOGD-treated cells. In the mouse model, the DHLF group with PVT1 overexpression had heavier brain tissue damage than the control group; apoptosis rate, reactive oxygen species level and caspase-3 enzyme activity were significantly higher in the lenti-PVT1 group than in the lenti-control group. Compared those in the lenti-control group, the total distance traveled, distance of action in the center, number of entering the center, average speed of walking reduced, and the distance of walking in the periphery and and peripheral walking distance increased in the lenti-PVT1 group. The dual luciferase reporter gene assay verified the possible binding sites between PVT1, miR-148a-3p and MKL1. In the animal model, cellular model, and blood samples of children experiencing DHLF, miR-148a-3p expression increased and MKL1 expression decreased. In experimental studies in vivo and in vitro, PVT1 and MKL1 expression increased, and miR-148a-3p expression decreased. Meanwhile MKL1 inhibitor CCG1423 reversed the apoptosis in neuronal cells. ConclusionPVT1 may adsorb miR-148a-3p to regulate the expression of MKL1, a downstream gene of miR-148a-3p, a mechanism promoting the apoptosis of neuronal cells in DHLF mice.

The emerging acute respiratory disease, COVID-19, caused by SARS-CoV-2 Coronavirus (SARS2 CoV) has spread fastly all over the word. As a member of RNA viruses, the glycosylation of envelope glycoprotein plays the crucial role in protein folding, evasing host immune system, invading host cell membrane, even affecting host preference. Therefore, detail glyco-related researches have been adopted in the Spike protein (S-protein) of SARS2 CoV from the bioinformatic perspective. Phylogenic analysis of S-protein sequences revealed the evolutionary relationship of N-glycosylation sites in different CoVs. Structural comparation of S-proteins indicated their similarity and distributions of N-glycosylation sites. Further potential sialic acid or galactose affinity domains have been described in the S-protein by docking analysis. Molecular dynamic simulation for the glycosylated complexus of S-protein-ACE2 implied that the complicate viral binding of receptor-binding domain may be influenced by peripheric N-glycans from own and adjacent monoers. These works will contribute to investigate the N-glycosylation in S-protein and explain the highly contagious of COVID-19.

The poly-ADP-ribose polymerase (PARP) superfamily consists of 17 members, which regulate many biological processes in physiological or pathological conditions, such as DNA damage repair, anti-viral responses, and development of adaptive immune cells. Among them, PARP14 is the biggest member, containing two RNA recognition motifs at the N-terminal, three macro-domains, one WWE domain, and one PARP domain at the C-terminal, which was reported to regulate IL4/STAT6 signaling in adaptive immune cells. However, whether PARP14 participates in regulating host inflammatory response remains unclear. In a previous study, we observed that virus infection and LPS treatment induced the transcription of Parp14. By comparing the primary macrophages derived from Parp14 KO and WT mice, we found that some inflammatory cytokines were significantly induced in KO macrophages. Still, the expression of Ifnb1 had no significant difference compared to the WT macrophages. RNA-seq analysis showed that the KO group had a more robust inflammatory response but a weaker innate immune response upon stimulation. We validated these results by performing a knockdown of Parp14 in RAW 264.7 cells. Moreover, the survival time of the KO mice was much shorter than that of the WT group upon LPS injection. Transcription factor enrichment analysis indicated that nuclear factor-kappaB1 (NF-{kappa}B1) may be the main reason for increasing the production of these inflammatory cytokines. As expected, the up-regulation was deleted upon the treatment of the inhibitor of NF-{kappa}B, JSH23. These data imply that PARP14 regulates inflammatory responses through the NF-{kappa}B pathway.